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Abstract 

The Earth's magnetic field effect on the azimuthal distribution of extensive 
air showers (HAS) of cosmic rays has been evaluated using a bulk of the 
Yakutsk array data. The uniform azimuthal distribution of the EAS event rate is 
rejected at the significance level 10"'"^. Amplitude of the first harmonics of 
observed distribution depends on zenith angle as A; ~ 0.2xsin and is almost 
independent of the primary energy; the phase coincides with the magnetic 
meridian. Basing upon the value of measured effect, the correction factor has 
been derived for the particle density depending on a geomagnetic parameter of 
a shower. 

The charged particle trajectories of the EAS developing in the atmosphere are 
distorted in the geomagnetic field. It results in stretched lateral distribution of particles in the 
plane perpendicular to the shower axis along the direction of the Lorentz force. The 
observation of this effect at the Yakutsk array in very inclined showers at zenith angles 9 > 
6(f has revealed the significant asymmetry [1] of the lateral distribution fiinction (LDF) of 
particles due to the deviation of muons proportional to the geomagnetic parameter, g, 
determined by the angle between the shower axis and the field vector: g = sin^^ / cos^0, 
where x = arccos(cos0 cos^h + sin0 sin^n cosa), 0h=14°, Ct - azimuth. 

The geomagnetic effect on the lateral distribution of particles was confirmed 
analyzing the showers in zenith angle range 20°<0 <60° [2]. The calculation results of the 
expected effect in the QGS model were used to analyze a shower of the highest energy 
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£'o>10 eV observed at the Yakutsk array [3]. In [4 a, b] it was shown that the showers 
initiated by extremely high energy photons have to show the north-south asymmetry. But the 
last effect has a threshold for the primary energy: Ey> 3- 10^^ eV. 

In this paper we draw attention to the fact that due to the LDF asymmetry depending 
on %, one can select the showers of equal energy and zenith angle in order to obtain the 
azimuthal modulation of the EAS event rate because of the different LDF stretch which 
causes the different observed primary cosmic ray energy and intensity. 

Fig.l shows the distribution of the number of EAS events above 5- 10^^ eV measured 
with the Yakutsk array in the period 1974-1995, zenith angles in the bins 20°-30°, 40°-50°, 
60°-70°. The number of EAS events in zenith angle intervals used in the analysis is given in 
Table 1. 



Table 1. The number of EAS events (n) 
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6 


7 
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e, deg. 


0-10 


10-20 


20-30 


30-40 


40-50 


50-60 


60-70 


70-80 


n 


25924 


69114 


83226 


66649 


37575 


16384 


6609 


1749 



The amplitude and phase of the first three harmonics of the distribution versus zenith 
angle are shown in Fig. 2, and versus the primary energy - in Fig. 3. The second and third 
harmonics amplitudes come up to the expected value for the uniform distribution of azimuths 
shown by a dashed line for the measured number of showers in each interval. Dotted lines 
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show the r.m.s. deviations from the expected amplitudes. The first harmonics amplitude 
essentially differs from zero in zenith angle intervals no. 3 to 7 given in Table 1. In these 
intervals one can discard the uniform distribution with error probability below 10'^'* basing 
on the probability for the uniform distribution to have the first harmonics amplitude Ai: 
P(>Ai)= exp(-nxAi /4). The first harmonics phase coincides with the magnetic meridian in 
Yakutsk. Almost vertical vector of the magnetic field in the array region (6h= 14°) results in 
the dominance of the first harmonics in the range 2(f <6 <7(f. For other arrays a situation 
may be different. For example, at the Tibet array where the field zenith angle is Oh = 45°, 
both the first and the second harmonics at 6 >50" are equally prevailing. At the Chacaltaya 
array ( 6h =88°j the second harmonics should dominate. 

To reveal the influence of the array's geometry on a modulation power, the azimuthal 
distribution of the event rate was used in a sample of showers registered by detectors 
arranged in a circle with radius 1.5 km around an array center (azimuth-symmetric part of the 
array). It turned out that the amplitude and phase of the first harmonics in the sample differ 
from the corresponding values of the general distribution less than the statistical errors. That 
is why we've assumed the azimuthal modulation independent of the array geometry. 

The data handling of the Yakutsk array used to be carried out fitting the parameters 
of the axially symmetric function, which approximates the charged particle density in 
detectors. The particle density at a distance 300 m from the shower core, paoo, found in this 
way, is connected to the EAS primary energy. Let's consider how paoo changes because of 
the geomagnetic scattering of particles with LDF of the form p, = c x f^, where r is a core 
distance in a perpendicular plane to the EAS axis. The density functions of positively and 
negatively charged particles are shifted with respect to a shower axis: 

p^^ = cx(((rxcos(p + df + r xsincp)''^ +((rXcos(p-df+r^xsin(p)''^^^)/2, 
where J is a displacement distance of the charges; (p a (0°-^360°). 
Averaging over a circle of r = 300 m we obtain.- 

(fsoo^ff^'^sooX (1 + (ri xd3oo/2f), 
where d«300; d3oo=d /300 m. On the average, the particle density in the detector increases 
in comparison with an expected density in the case of absence of the geomagnetic field (i.e. g 
=0). For showers arriving from the north the observed particle densities in detectors are 
higher than in "southern" showers of the equal energy with the equal zenith angle. 

When we select showers with the same observed densities pjoo, the primary energy 
of showers arriving from the north is less. It gives the decrease in a corresponding event rate 
as it is seen from Fig.l, because of the number of the registered showers at the Yakutsk array 
diminishes as the energy decreases in the region 10^^<£'o<10^^ eV. The amplitude of the first 
harmonics Ay is defined by a relative value of density change measured for the southern and 
northern showers: 2xAi = AJ / J = yxApsoo/ psoo, where J is the EAS event rate; 7= InJ / In 
psoo- The zenith angle dependence of A7 is well described by Q.lxsin (Fig.2), but the 
variation of y with zenith angle leads to the ratio Apsoo/psoo, which is not described by the 
simple dependence. 

It is known that the zenith angle boundary, separating showers where electron-photon 
or muon component dominates, is about 50" at the sea level, r = 300 m [5]. Fig.2 shows that 
the geomagnetic effect in EAS is observed in both these regions. A distance of the charge 
shift determines the LDF asymmetry caused by the geomagnetic effect. The ratio ^ of the 
largest distance to the least one corresponding to the fixed particle density p^'^^soo is: 
^ ~ \+(\ + ri/2)x 300 ~ l+Apsoo/psoo for r] = 3. So, the measurements of the LDF "oval" can 
be used to estimate the change of psoo due to the geomagnetic effect. As it was shown in [1], 
the linear coefficient between the LDF asymmetry and the geomagnetic parameter is d^dg= 
0.1 ± 0.04 for very inclined showers. Using this value for all zenith angles we obtain the 
correction factor to the density psoo (Table 2). 
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Table 2. Correction factors to the measured density pjoo 







20 


30 


40 


50 


60 


"300/(^300 


Northern 
showers 


0.94 


0.92 


0.88 


0.82 


0.72 


Southern 
showers 


0.99 


0.96 


0.93 


0.87 


0.78 



As it is seen from Fig.3 the first harmonics amplitude is almost independent of Eo in 
the energy range where there is sufficient number of events, so the factor dependence on 
energy can be neglected. The recalculation from psoo to Eo at the Yakutsk array is carried out 
using the method of connecting P300 in inclined and vertical showers through the equi- 
intensity lines in spectra at the various zenith angles. Because of these intensities correspond 
to the same energy of the EAS primaries and the geomagnetic corrections can be omitted in 
vertical showers, then the procedure of Eo estimation correctly takes into account the zenith- 
angle dependence, but averages the variation of energy with azimuth. It leads to a small 
systematic error proportional to in estimation of the primary energy. The attenuation 
length of p3oo obtained by equi- intensity lines is overestimated by ~10% in comparison to a 
value for g=0. 

When the geomagnetic corrections are taken, the EAS size spectrum index for psoo 
(P600 for showers registered with detectors spaced by 1 km distance) increases due to the 
augmented with P300 fraction of the inclined showers detected by the array. The resultant 
corrections to the size spectrum parameters are given in Table 3. 

T able 3. The spectrum index (y) and intensity (J) co rrections 



Eo,eV 


10^^ 


3x10^^ 


10^^ 


Ay 


0.01 


0.07 


0.24 


AJ/J 


0.16 


0.19 


0.28 
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Captions to Figures 

Figure 1. Relative numbers of EAS events in zenith angle intervals (shown at the right) for 
showers above 5-10^^ eV. The function l+AiXcos( a -Ui) is shown by dashed lines. 

Figure 2. The amplitudes of the first three harmonics (upper panel) and the phase of the first 
harmonics (lower panel) versus zenith angle. Statistical errors are shown by the vertical bars, 
horizontal bars show the angular bins. Symbols are: triangles - the first harmonics; circles - 
the second harmonics; squares - the third harmonics. The approximation O.lxsin 6 is shown 
by the dash-and-dot line. The expected amplitude of harmonics for the uniform distribution 
in azimuth is given by the dashed Une. The standard deviation of expected amplitudes is 
shown by the dotted Unes. 

Figure 3. The amplitudes of the first three harmonics and a phase of the first harmonics as a 
function of energy. The symbols are the same as in Fig.2. 
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